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I. INTRODUCTION 



A simple relativistic version of the Schrodinger equation is the Spinless Salpeter Equa- 
tion (SSE). For the one-dimensional case we have 

^-dl + m2 ^(x) = (E- V{x)) ^(x), (1) 

where m is the mass of the particle, V{x) is the potential interaction, E the eigenenergy 
of the stationary state d^. = = —j? and p is the relative momentum of the 

particle = c = 1). p and x are conjugate variables. The differential operator of the 
Schrodinger equation is well defined because it is a second derivative. To solve a physical 
problem, we must just solve an ordinary eigenvalue differential equation. The situation is 
more complicated with the SSE because the associated differential operator is a nonlocal 
one. Its action cannot be calculated directly from its operator form. Indeed, its action on a 
function /(x) is known only if j{x) is an eigenfunction of the operator d^. In this case we 
obtain 

^-dl + m2 f[x) = V-a + w? fix), (2) 

where a is the corresponding eigenvalue of d"^. That is why we need first to rewrite the 
SSE into a form easier to handle. Since the operator is a nonlocal one, this form could 
be an integral equation. This have been done for the three-dimensional case in Refs. [|T],0|. 
We present the one-dimensional corresponding form in the next section. With the method 
used to obtain this form, it is possible to rewrite the SSE as an integro-differential equation 
(see Ref. for the three-dimensional case). But the kernel is really complicated and the 
resulting equation seems to be very difficult to treat. We will use, here, another method to 
obtain the solution of the equation. 

To solve the relativistic Coulomb problem we do not solve any differential equation. 
We calculate the action of the square- root operator on the functions x"e~^^. Because the 
result is analytical and because the wave functions of the Coulomb problem with a hard 
core interaction are an exponential multiplied by a polynomial, which is also the form of 
the Schrodinger and Klein-Gordon wave functions, a complete solution of Eq. with 
V{x) (xl/x and x > 0, can be found. 

The paper is organized as follow. In Sec. ^ we give some useful mathematical results 
concerning the square-root operator. In Sec. [ITTt we solve the one-dimensional Coulomb 
problem with a hard core interaction. In Sec. |rV|, we compare our results to those obtained 
with the Schrodinger equation and with the Klein-Gordon equation [^|[. At last, we 
give our conclusion in Sec. [V]. 

II. MATHEMATICAL FRAMEWORK 

In this section we give some results concerning the square-root operator which we use to 
solve the Coulomb problem supplemented by a hard core interaction. 
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A. Integral representation of the square-root operator 

To obtain the integral representation of the square-root operator we use the Fourier 
transform of the one- dimensional delta function. We have 

J-dl + m? m{x) = — dpdq Jv^ + rn^ e'^^^"^)*' *(g). (3) 

Extracting the operator — + w? and integrating over the momentum p (see Ref. we 
obtain, 

I 1 / \ /■+°° 

* TT ^ ^ J-oo 

= -(-dl + m^) dq Ko{mq)[^{x + q) + '<iJ{x - q)], (4) 

TT ^ ^ JO 

where Kq{x) is the modified Bessel function of order p. 952]. 



B. Invariant space functions of the one-dimensional square-root operator 

In this section we calculate the action of the square-root operator on the functions x"e~^^. 
We obtain that the result is equal to a polynomial of order n, Mn{m, P, x), multiplied by the 
same exponential. Thus the space of functions Pn{x)e~^^ is the invariant space functions of 
this operator. Using formula (Hf) we have 



-dl + m? a;"e-^" = - i-dl + m") e"'^" / dq Ko{mq) {x + g)"e-'^« + {x - qfe'^'' 

TT ^ ^ JO '- 



(5) 



This leads to |, p. 712] 



The coefficients Gk{Tn,f3) are given by 

= fl^S ( (^^TTT)^^ («= + U/2; * + 3/2; + i-)%3 ^ . 

(7) 

where F{a, (3;'y; x) is the hypergeometric function p. 1039]. Performing the derivation 
in Eq. (|^) and rearranging the obtained relation, we have: 



-dl + m2 a;"e-^" = ^ e''^" E ( ^ ) | - Gk{m, j3) + 2/3A;G'fc_i(m, /?) 



^ fc=0 



A;(A;-1)G,„2K/?)| x"-^ (8) 
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It is possible to write the coefficients Gfc(m, /5) into a more useful form. This form will 
allow us to find a recursion relation between the coefficients Gfc(m, /?) and to simplify the 
expression (P). We will be able to construct the polynomial, M„(m, for each value of 
n. We have the relation [jlO|, p. 562] 



(l-2a)/4 /I - T 

F(a, 1/2; a + 1/2; -x) = r{a + 1/2) _ P'{%" ' 



X 



1 + x 



(9) 



with X > 0. The functions Py{x) are the associated Legendre functions for x real and 
\x\ < 1. With this relation, we find that 



T{kf 



-(3/m) 



^ (^2_^2)(2A:-l)/4 

Now, using the recursion relation of the associated Legendre functions 0, p. 1005] 



-2(/x+l) 



X 



\/l — x- 



,Pl^+\x) + {ii-u){ii + u + l)PI^{x) 



(10) 



(11) 



and the explicit expression of P^i/2 {x) and P_^/2{^) a-> P- 1008], we can write the following 
relations 



Gk+2{m,(3) 



{k + If Gk{m, (3) - {2k + 3)/3 Gfc+i(m, (3) 



with 



Go{m,P) 
G,{m,(3) 



TT 



(m2 - /32) 



(12) 

(13) 
(14) 



Al last, one can find, using Eq. (P^), that the general coefficient of the sum of the Eq. 
becomes 



n 
k 



(3Gk-i{m,P)-{k-l) Gk-2{m,P) 



with k > 1. 



(15) 



And with this form, a recursion relation for F/; „(m, (3) can be easily found. Thus to conclude 
this section, we are able now to rewrite Eq. (H) into a simple form: 



Fo,„(m,/?) = v'm2-/52, 
n/5 



^Ffc,„(m,/?) X 

lk=0 



n—k 



-I3x 



Fi^n{'m,(3) 
Fk+2,n{m,(3) 



Vm2-/32' 

n — A; — 1 
(A; + 2) (m2 - /?2) 

(2A;+l)/5F,,+i,„(m,/3) 
n + l 



k,n+l 



n + 1 — k 



Fk.n- 



(16) 

(17) 
(18) 

(19) 
(20) 



We can see that we obtain the expected relation (from Eq. (^) for n = 0. And thus we 
see that we must have [3 < m. With the relations (|T7|) -(pO|) the polynomial Mn{m, P,x) is 
completely defined and we can construct it for each value of n. This result will allow us to 
find, with few calculations, the solution of the one- dimensional relativistic Coulomb problem 
with a hard core interaction. We give below the polynomials, as an example, for n = —>■ 4. 



Mo{m,P,x) = S, (21) 



Mi{m,/3,x) = S x + ^, (22) 



2/3 


X - 


2 

m 


S 




3(3 


x' 




'J 


S3 


A(3 


x^ 




S 


S3 



with 



(23) 

^+^, (24) 
M4(m,/?,x) = Sx' + ^x'-^-^x' + 1^ ^ - ^ (m^ + 4/?^) , (25) 



S = y^m2 - /32. (26) 



Note that these last relations can be simply checked by acting the square-root operator on 
each side of Eq. (|TB|). For n = 1, we see that we have an identity if we use the relation 
for n = 0. Now knowing these two relations we see that the relation for n = 2 is also an 
identity, and so on for each value of n. 



III. THE ONE-DIMENSIONAL RELATIVISTIC COULOMB PROBLEM WITH A 

HARD CORE INTERACTION 

The equation to solve is 

^-dl + m? ^(x) = + ^) ^(x). (27) 

We just consider here the case x > (we will discuss after the extension to the whole x 
axis). Physically this means that we have a hard core interaction for x < 0. Then the 
wave functions will possess the following asymptotic behavior: \E'(x) = for x < and for 
X = +00. Suppose that the wave functions have the following form 

n 

\E'(x) oc ^ 7fc,„ x^ e"^^ for x > and n = 1, 2, 
fc=i 

^(x) = for X < 0. (28) 

We do not consider the normalization of the functions here. Thus replacing Eq. (^) into 
Eq. ( ^71) and using Eq. ([I6|) we obtain 

n k n n 

E 7fe,n E FpA^, P) ^^"^ = ^ E 7fc,n x^ ^kY, 7fc,„ x^-^ (29) 

k=\ p=0 fc=l k=\ 
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Now equaling order by order we will determine the solution. The term of order n gives: 



E = Fo,n{m,f3) = ^m^-(3^. (30) 
From the term of order n — 1, we have: 

K = Fi,4m,f3), (31) 

which leads to, 

^ (32) 

ijy'l + (K/n)2 

We can remark that we have as well the necessary relation (3 < m. We are now already able 
to determine the energy spectrum. Using Eq. (|30D and Eq. (B^) we have: 



E = . ^ . (33) 

To obtain a complete solution, we must now find all the 'yk,n, and prove that the system of 
equations which gives these quantities has always a solution. Obviously we can fix 'jn,n = 1- 
We see that the term of order n—j determines the coefficient 'yn-j+i,n if the previous 7fc,„ are 
known. Beginning with the term of order n — 2, we obtain directly 7n-i,n- And now we can 
get 7„-2,ra from the term of order n — 3. The independent term will fix the last factor 7i^„. 
Thus we have a triangular system of n — 1 algebraic equations with n — 1 unknowns. This 
system will possess always a solution if the determinant of the coefficient matrix is non null. 
As this is a triangular matrix, the determinant is the product of the diagonal elements. The 
expression of these elements is k — Fi „_j(m, /5) which is equal to jP/S. These quantities 
are always non null since j > 0. The general form of 7fc^„ is obtained from the term of order 
n — j — 1. We have 

S •'^^ 

1n-j,n = —3^1 ln-k,n Fj.k+l,n-k{m, (3). (34) 

k=o 

We can inverse the summation to finally obtain 

j-i pi-i Pj-i-i 

7n-i,n = X! I] ••■ I] F{n,Pu3)F{n,p2,pi)...F{n,pj,pj^i), (35) 

with 

S 

F{n,kJ) = — Fj^k+i,n-k{m,(3). (36) 

For the summation in Eq. (^), we must use the following rule: If in a summation over p^, 
a being arbitrary, the bound Pa-i — 1 is negative, all the F{n,k,j) containing the indices 
Pi3>a are equal to 1. With the formula (0), we are able to construct the wave functions for 
the Coulomb problem with a hard core interaction. As an example we give the three first 
one: 
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^(x) oc X Qn{m, K, x) e"^"", (37) 

with P given by Eq. (0), and 

Qi{m, K,x) = 1, (38) 

2 

771 

Q2{m,K,x) = X - (39) 
Qsim, K,x)=x^-^ x + [(3^ + m2) , (40) 



with S defined by Eq. (p6|). Again, we can perform a simple verification by putting these 
solutions into Eq. (^) and using Eq. (plGl). 

Contrary to the Schrodinger or Klein-Gordon equation, the extension of the solution 
to the whole x axis is really more complicated. We can try to use exp(— /3|x|) instead of 
exp(— /3x) in our solution. But the situation is quite more difficult. Indeed, the construction 
of the solution was based on the fact that exp(— /3x) was an eigenfunction of the square-root 
operator and that the invariant space functions of this operator was Pn{x) exp(— /5x), where 
P„(x) is a polynomial of order n. But it is easy to show, with Eq. (H), that 

I 2772 

J—dl + m'^ exp(— ?77|x|) = — Ko{m\x\). (41) 

This is non null as this is the case in the Eq. ([I6|). Thus exp{—(3\x\) is not an eigenfunction 
of the square-root operator and Pn{x) exp{—(3\x\) is not an invariant space function of this 
operator. So it seems that the pure Coulomb problem has quite different solutions for the 
wave functions and certainly for the spectrum. 



IV. DISCUSSION 

The one-dimensional Coulomb problem has been treated by many authors, both nonrel- 
ativisticly and relativisticly f^,]^. But in these works the whole x axis is considered. 
As a consequence, the ground state gives some difficulties. 

In the nonrelativistic case the solution is: 

\l/(a;) = X exp(— 7cr77|x|/77) L^}^i{2Km\x\/n), (42) 

E = m with 77 = 1,2,... (43) 

But we see that for n = 1, the wave function have a node at the origin. So this is not the 
wave function for the ground state. In fact it is found to be infinitely bounded and the wave 
function is a delta function 

In the Klein-Gordon case the solution is: 

^{x) = x' exp{-(3\x\/2) l2,(/3|x|), (44) 
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with, 

(3 = 2mK/^{n - I + Sf + with n=l,2, (46) 

and, 

S = \{1 + i) = \{i±^/i^a;^). (47) 

Thus we see that we have two distinct solutions according the sign for S. Actually, spectrum 
with the minus sign for S is not acceptable. Indeed, a reason is that, for n = 1, when we 
perform the limit k — 0, we obtain = 0. This mean that the particle is still bounded 
when the interaction vanishes. Thus the problem for the ground state persists (see Ref. Q 
for a complete discussion). 

In this paper we do not consider the whole x axis and we have no problem with the 
ground state. We consider a hard core interaction, for x < 0, which gives \E'(x) = in this 
region. Thus xexp(— /5x) is the wave function for the ground state. Actually the purpose of 
these work was to solve a particular kind of differential equations with a difficult to handle 
nonlocal operator. Indeed, any analytical solutions are known for the spinless Salpeter 
equation. Thus we do not discuss the problem of the ground state of the one-dimensional 



Coulomb problem. In Sec. [ITI| , we have shown that the extension to the whole x axis is not 
easy. Moreover, the ground state problem could persist. 

To compare our result to the results of previous works, we can consider the Schrodinger 
and the Klein-Gordon equation for the Coulomb potential supplemented by a hard core 
interaction. The spectra and the wave functions remain unchanged but the ground state 
problem has disappeared. In the three cases we have the same kind of wave functions: an ex- 
ponential (with different arguments) multiplied by a polynomial (with different coefficients). 
For the spectrum we have, in the limit of small k: 

^sch = (l - 1^) ' (48) 

-.—1 ( ^ ^ ^ ^Kj \ . ,„v 

^KC = -„^l-^--, + -j. (49) 

i Sun , , 

^^--"(1-^ + ^). (50) 

Thus we see that in the expansion of the Salpeter spectrum the term in v? is missing 
compared to the Klein-Gordon spectrum. So the difference between these two spectra is 
rather important. For an electron in an electromagnetic Coulomb potential, the splitting is 
about 10-3 eV. 

Another characteristic of the spinless Salpeter spectrum is that k can grow up without 
limit. This could come from the fact that we have another kinetic operator than in the Klein- 
Gordon equation and that the result could be quite different. But the main explanation is 
certainly that we do not solve the real Coulomb problem and that this spectrum could be 
different contrary to the Klein-Gordon equation, which keeps the same spectrum in both 



cases. Indeed, for the SSE, there exists a limit value for k in three dimensions 11 . 
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V. CONCLUSION 



The purpose of this work was to find an analytical solution of a particularly kind of 
differential equations containing a nonlocal differential operator. The equation considered 
in this paper, the one-dimensional spinless Salpeter equation (SSE), is a simple relativistic 
version of the one-dimensional Schrodinger equation. The SSE is not a marginal equation. 



For three dimensions, this equation comes from the Bethe-Salpeter equation [12,13, p. 297], 
which gives the correct description of bound states of two particles. Moreover, despite the 
presence of a so particular operator, the SSE is often used in the potential models (see for 
instance Refs. []T^|-Q), which give a phenomenological description of hadrons. 

To find this analytical solution, we calculate, in Sec. |II B| , the action of the square-root 
operator on a polynomial multiplied by an exponential and we show that this constitute the 
invariant space functions of this operator. To be able to perform this calculation, we have 
constructed, in Sec. 11 A, an integral representation of the square- root operator. In Sec. |T|, 



we have obtained, without solving any differential equation, a complete solution of the SSE 
with a Coulomb potential and a hard core interaction. This last interaction is introduced to 
keep the particle in the x positive region. We remark that the SSE wave functions have the 
same form than the Schrodinger and the Klein-Gordon wave functions. We remark also that 
the splitting between the SSE and the Klein-Gordon spectrum is rather important. Indeed, 
it is of the same order of the first relativistic correction given by these two equations. 
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